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I Recent data from Relativistic Heavy Ion Collider (RHIC) on rapidity and 



transverse momentum distributions of hadrons produced in ultra-relativistic 
reaction of both Au+Au and p+p are compared to predictions of the HIJING 



> 

I and RQMD models. The original default mini-jet scale po = 2 GeV/c and 
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, the multiplicity with energy. RQMD model without mini-jet leads on the 

other hand to a too slowly increasing multiplicity with energy. Therefore, we 

O I study what variations of po and dE/dx are required in HIJING to account 
^3 ' 



for the observed Npart and ^/s dependence of the global dNch/dy and dEr/dy 
observables as well as the jet quenching pattern out to ~ 8 GeV. We show 



X 

\ that a slight increase of po from the default 2.0 GeV/c at ^/s = 130A GeV to 

2.2 GeV/c at ^/s = 200 A GeV is sufficient to account for the bulk observables. 
Jet quenching of vr*^ above px > ^ GeV/c at 200A GeV is found to be well 
reproduced in central Au + Au by the original HIJING default prediction. 
However, the effective surface emission in the HIJING model formulation of 
energy loss over-predicts the quenching in more peripheral collisions. Neither 
HIJING nor HIJING/BB can account for the observed anomalous excess of 
moderate pT < ^ GeV/c baryons. The agreement of RQMD with data in 
this pt region is shown to be fortuitous because, without mini-jets, it fails to 
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reproduce the p + p spectrum. 

PACS numbers: 25.75.-q; 25.75.Dw; 25.75.Ld 



I. INTRODUCTION 

Qualitatively new phenomena have been discovered [1,2] in ultra- relativistic nuclear col- 
hsions at total center of mass energy Ecm = ISOA — 200A GeV at the Relativistic Heavy 
Ion Collider (RHIC) and interpreted as evidence for the formation of ultra-dense QCD 
matter. The extensive systematics of both bulk global observables and high transverse mo- 
mentum phenomena were reported by PHENIX [3-13], STAR [14-23], PHOBOS [24-31] 
and BRAHMS [32-36] . RHIC data are now available not only for Au + Au collisions but 
also for p -|- p with high statistics, to accurately calibrate the magnitude of nuclear effects, 
and also d-\- Au [13,23,31,36] needed as a control experiment to separate initial versus final 
state interaction effects. 

In this paper we compare predictions of the HIJING [37] and RQMD [38] models to 
data. These nuclear collision event generators were developed a decade ago and continue to 
be useful tools for detector design and interpretation of experimental results because they 
simulate complete exclusive event characteristics at SPS and RHIC. HIJING, incorporating 
pQCD mini-jet production, is constrained to reproduce essential features oip-\-p data over 
a wide energy range. RQMD is designed to simulate only soft multi-particle production 
at lower energies, but is of interest since it incorporates a model of final state interaction. 
No attempt is made here to review the many other models developed since the release of 
HIJING and RQMD. We refer the reader to Refs. [39-76] for a broader perspective. 

Our goal is to test specific predictions within HIJING and RQMD models, made well 
in advance of the data. The predictions involved estimating key physical parameters con- 
trolling the dynamics. An important aim of this paper is to investigate what adjustments 
of those parameters may be required in light of the new data. The two main parameters 
of HIJING that we concentrate on in this paper are (1) the separation scale, po between 
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the perturbative (pQCD) mini-jet processes and the phenomenological soft (beam jet frag- 
mentation) processes and (2) the energy loss, dE/dx, of high transverse momentum partons 
propagating through the dense (quark-gluon plasma) medium produced in the reaction. 

The mini-jet separation scale, po ~ 2 GeV/c, assumed in HIJING [37] controls the 
i/i dependence of the bulk multiplicity and transverse energy observables as well as their 
centrality dependence on the number of wounded nucleon participants N^art- RQMD [38] 
and its UrQMD [39,71] extensions assume in effect that pq ^ oo and therefore neglect the 
power law tails due to pQCD mini-jets. The default constant value of po 2 GeV/c in 
HIJING was fixed by fitting pp data through Tevatron energies [37], and it was assumed 
to be independent of both ^Js and A in order to predict multiparticle observables in AA 
that extrapolate down accurately to ^4 = 1 to reproduce experimentally known pp data. In 
RQMD, there was no attempt to fit collider energy pp data. Our approach differs in this 
major respect from many recent models that fail to account for multi particle phenomena 
in p -|- p collisions. Our philosophy is that due to the myriad of dynamical complexities 
already displayed in "elementary" p + p reactions, any model proposed to explain A + A 
collisions must extrapolate accurately down to ^4 = 1. This is because one of the very few 
experimental knobs in heavy ion reactions is the variation of the impact parameter through 
the participant number dependence of observables. In peripheral collisions the observables 
necessarily approach their value in p -|- p collisions. 

Recently, several models were developed that challenge the assumption that the separa- 
tion scale po is independent of both ^/s and A as in HIJING, Refs. [57,58] generalized the 
mini-jet scale by allowing it to vary dynamically by introducing the hypothesis of final state 
saturation (FSS) of the produced mini-jet density per unit area. That saturation scale was 
predicted to be Ps{s,A) = 1.1 GeV {,/^/200f^^^{A/200)^-^^\ This hypothesis leads to a 
specific prediction for the initial parton distribution in A + B collisions. With the additional 
hypothesis that this initial condition evolves according to local thermal hydrodynamics, the 
Pt integrated global observables in central Au + Au were shown to be well account for. How- 
ever, the breakdown of hydrodynamics and the smallness of p^ (p^ < 1 GeV/c) in peripheral 
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collisions prevents the model from predicting correctly peripheral A + A collisions and p+p. 

An alternate initial state saturation (ISS) hypothesis for the variation of the separation 
scale was introduced in Refs. [59-63,77,78]. In this picture, po is replaced by a gluon satura- 
tion scale Qs{s,A) ^ 1.4 GeV (v/s/200)°-^(A/200)V6. This saturation scale is considered to 
be the boundary of the classical Yang-Mills field domain. Instead of hydrodynamics, local 
parton hadron duality is assumed to predict low pt ~ Qs integrated bulk global observables. 
The normahzation of Qs was fixed by fitting the observed central 130 A GeV Au-|-Au rapidity 
density. ISS was found to be more successful in describing the lower participant number and 
^/s dependence of the rapidity density at RHIC than the FSS model due to a particular low 
(5^ ~ 1 — 2 GeV^ dependence of the gluon structure function and fine structure coupling. 
While prehminary extensions [62] of this ISS model to the pt > Qs regime could fit central 
130A GeV Au + Au without final state interactions, the most recent data on d + Au reac- 
tions [13,23,31,36] rule out a particular extension [63] of the ISS model at mid rapidity in 
the range 2 < < 10 GeV/c. 

While neither ISS or FSS saturation models can describe simultaneously the global low 
Pt and the hard high pt observables in both p-\-p and A-\- A collisions, they both provide 
strong motivation to test the effects of variations of po with both s and A in HIJING. 
With this motivation, the first part of this paper will be to investigate, within the HIJING 
model, the effect of relaxing the {s,A) independence of the default constant po = 2 GeV/c 
assumption. Importantly, we seek to do this while retaining our philosophy that critical 
features of multiparticle production in p + p should be accounted for simultaneously in the 
same model. Our first conclusion will be that the global 130-200A GeV Au -\- Au data can 
be well accounted for by allowing a rather modest 10% (A independent) enhancement of po 
from 2 GeV/c at 130A GcV to about 2.2 GeV/c at 200A GeV. 

The second part of this paper focuses on jet quenching [79]. Jet quenching is one of the 
most striking new phenomena [6,11,17,20,22,30] discovered at RHIC. This effect was not 
observed previously at lower energies. In fact, at SPS a strong (Cronin) enhancement of 
high Pt 7r° was observed [80] . 
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The first quantitative predictions of jet quenching with HIJING [37,81] assumed a default 
constant dE/dx = 2 (1) GeV/fm gluon (quark) energy loss and was implemented by a 
simple string flip algorithm assuming a mean free path of A = 1 fm. As recently reviewed in 
[40,82], there has been considerable progress since that time in computing the dependence 
of medium induced radiative energy loss, AE(pt, L, pq), in QCD [45-47,52,53,82]. The 
observable consequences of the dependence of AE on the jet energy, its propagation length, 
and the evolving parton density have been explored in pQCD models in which the plasma is 
assumed rather than calculated dynamically. While the jet quenching algorithm in HIJING 
is much more schematic, the model provides a useful theoretical laboratory to study its 
observable consequences in the dynamical medium that it creates. In the second part of 
this paper, we therefore explore what modifications of the default HIJING assumptions are 
required in light of the new data. 

As we show below, the default constant energy loss in HIJING accounts remarkably well 
for the high px 7r° suppression pattern in central Au + Au at 200A GeV. In fact, HIJING 
also accounts [83] for the enhancement observed [80] at SPS that is otherwise puzzling 
according to pQCD estimates [84] . However, we show below that HIJING fails to account 
for the observed centrality dependence as well as for the anomalous baryon enhancement 
observed up to 5 GeV. HIJING/B [85-87], with its implementation of baryon junctions 
[88,89], was tested to see if this mechanism could account for the baryon anomaly (see 
also [48]). However, the present version (HIJING/BBvl.lO) fails to account for the large 
transverse slopes of anti-baryons and does not reproduce the "baryon lump" at moderate 
Pt in the nuclear modification factor. 

While RQMD [38] does not contain mini-jets, we investigate its predictions because it is 
one of first models to simulate final state transport dynamics of pre-hadrons and hadrons 
without assuming local equilibrium as in hydrodynamics. It was able to reproduce the 
directed and elliptic collective flow systematics observed in Pb+Pb at SPS (17A GeV). 
However, as we show below, the absence of hard pQCD processes leads to a much too weak 
beam energy dependence at RHIC energies. It therefore also fails to account for the power 
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law tails of the p+p spectra at 200 GeV. We refer to Refs. [71,72] for a review of application 
of RQMD and URQMD applied to reactions at AGS and SPS energies. 

While not addressed directly in this paper, we call attention to the recent AMPT trans- 
port model [65] that incorporates mini-jet production and extends HIJING by including both 
parton cascading and hadronic final state interactions. AMPT is under extensive develop- 
ment and has been tested on a number of important RHIC observables [65-68] . However, 
problems with covariance of numerical solutions involving ultra-relativistic parton cascad- 
ing [69,70] require very high parton subdivision techniques which are unfortunately beyond 
present computer power to solve with AMPT. 

In order to bypass current technical difficulties of predicting bulk collective phenomena 
via transport theory, hydrodynamic models have also been extensively applied [41,90-97]. 
The central simplifying dynamical assumption is that perfect local equilibrium is established 
and maintained throughout the reaction. Therefore non-viscous hydrodynamics together 
with a Cooper- Prye statistical freeze-out prescription [98] are used to compute the expansion, 
hadronization, and subsequent expansion until freeze-out. No attempt is made in such 
models to compute the initial condition, but rather the initial entropy and baryon density 
are fit to the measured rapidity distributions. While such models cannot predict beam 
energy dependence of observables, they do predict striking collective flow phenomena and 
their dependence on the QCD equation of state (for a recent review of hydrodynamics at 
RHIC see Ref. [41]). The first attempt at a hybrid combination hydrodynamics and jet 
quenching was proposed in [49]. Recently an important step forward is the development of 
a consistent 3-l-lD hydrodynamical approach including QCD jet quenching [97]. 

Unlike hydrodynamics [41,97] or parton transport models [66,69], neither HIJING nor 
RQMD can predict the large amplitude elliptic flow observed [15,21,7,29] at RHIC. Elliptic 
flow is especially sensitive to early partonic final state interactions beyond the capability of 
these models. However, azimuthally integrated inclusive spectra are still interesting and can 
be addressed by models considered here. Wc use the following versions of HIJING vl.37 
[37], HIJING/BB vl.lO [85,86] and RQMD v2.4 [38] for the computations reported below. 
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II. CHARGED PARTICLES DISTRIBUTIONS AND THE MINI-JET SCALE 

Recent measurements of the rapidity density of charged particles in Au+Au coUisions over 
the range of total nucleon-nucleon center of mass (cm.) energy ^/sNN=56 GeV - 200 GeV, 
have been reported [4-6,8,16,17,26-28,33]. Within the errors, an approximatively logarith- 
mic rise of charged particle rapidity density per participating baryon pair {dNch/dr)/0.5Npart) 
with y/sNN is observed over the full range of coUision energies [28]. In Refs. [99,100] the cen- 
trality dependence of this observable was proposed as a test of the nuclear enhancement of 
the mini-jet component as well as whether gluon saturation is reached at RHIC energies. The 
predictions of different models varied prior to the data greatly from dNch/dy ~ 700 — 1500 
at mid-rapidity for Au-|-Au central collisions [39]. 

The predictions of HIJING (Figs la,b) and HIJING/BB (Figs lc,d) with (y) or without 
(n) effects of quenching (q) or/and shadowing (s) are presented in Fig. 1. The data from 
PHOBOS [26-28] and BRAHMS [33] experiments at ^^^^^=130 GeV and ^5^^=200 GeV 
are shown for comparison. In all cases, both quenching and possible parton shadowing in- 
fluence the predicted pseudo-rapidity distribution, dNch/drj, in this energy range. However, 
these effects work in opposite direction and thus partially cancel each other. Without shad- 
owing as assumed in the default HIJING model, the flux of mini-jets with pt > 2 GeV/c 
is too high and dNch/drj is overestimated. Even with the larger shadowing at the smaller 
X at Y^s/vjv=200 GeV, the gluon density enhancement resulting from mini-jet energy loss 
leads to a 10-20% overestimate of the charged particle rapidity density in going from 130A 
to 200A GeV in Fig. lb. A similar tendency is seen in the BB version of HIJING which, 
however, better accounts for the width of the rapidity distributions. The width is sensitive 
to the nuclear fragmentation physics, especially baryon number transport from the beam 
rapidities. HIJING/BB can better account for nuclear fragmentation by introducing the 
greater baryon stopping power through the baryon junction mechanism. 

These data are also consistent with the initial state saturation ISS model [33,59,60]. 
However, the EKRT flnal state saturation model [57] tends to over-predict the width of the 
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rapidity distribution. 

The energy dependence of the particle multiphcity is more easily seen in Fig. 2 where 
the central rapidity density per participant pair vs ^/s of both HIJING and RQMD are 
compared to data. The PHOBOS data for the central (0-6%) Au+Au coUisions are from 
[26-28]. The data for pp and pp are from [101-104]. Beginning about ^/sNN=100 GeV, the 
central Au+Au collisions show a significantly larger particle density per participant pair than 
in inelastic pp coUisions. The energy dependence predicted by HIJING is strikingly different 
than that predicted by RQMD. While RQMD predicts a very small increase over the range 
Y^sjv7v=56 GeV - 200 GeV, HIJING predicts an increase of more than a factor of 1.3, which 
continue up to the highest energy calculated. This increase in HIJING is due to copious 
mini-jets production in A + A collisions. RQMD fails to describe the trend of data because 
it misses the rise in multiplicity due to mini-jets. The predictions of both HIJING and 
HIJING /BB models are in better agreement with the data when the effects of both quenching 
and shadowing are included. Note that with default energy loss (dE/dx=2 GeV/fm) and 
constant po=2 GeV/c, the energy dependence obtained with HIJING is too rapid and the 
curves in Fig. 2 are for a reduced effective energy loss dE/dx = 0.5 GeV/fm. One observes, 
however, that even this smaller energy loss still leads to a more rapid dependence on energy 
than seen in the data. The RQMD v2.4 curves were obtained using the cascade mode and 
taken into account its rescattering and color ropes effects with their default parameters. 

Another global probe of the dynamics is the transverse energy per charged particle, 
dET/df] [4]. This distribution is sensitive to PdV work done by the plasma in hydrodynam- 
ical models. In HIJING its value depends again on the assumed shadowing and energy loss 
as seen in Fig. 3. In part (a) the total charged plus neutral transverse energy distribution is 
shown. In part (b) the contribution from only charged particle is shown. The results are for 
central (0-5 %) Au-|-Au coUisions at ^/s^—130 GeV. Both versions HIJING (yq,ys) and 
HIJING/BB seem to account better for the observed dEr/drj ^ 500 GeV than RQMD. This 
is again due to the absence of mini-jets in RQMD model. 

We recall that FSS saturation model tend in contrast to over-predict [57] by a factor 
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2-3 the transverse energy because the saturation scale Ps ~ 1 GeV/c is significantly smaller 
than the default po = 2 GeV/c needed in HIJING to fit p data. Therefore, FSS requires 
reduction of the initial transverse energy due to longitudinal hydrodynamic work. The same 
general tendency of over-predicting the transverse energy is found in classical Yang Mills 
simulation of A-l-A [43] . However, no detailed predictions of the transverse energy have been 
made within the KLN version [59-61] of ISS models. 

Another important difference between the predictions of models is in the rapidity de- 
pendence of the transverse energy per particle (see Fig. 3c and 3d) . While RQMD predicts 
a relatively constant value between —2.5 < f] < 2.5, both the numerator and denominator 
disagree with the data. HIJING gives on the other hand a rather strongly peaked distri- 
bution at mid-rapidity. This peaked distribution in HIJING is due to the locahzation of 
mini-jet production to central rapidities. Hydrodynamic models [41] generally assume a 
uniform boost invariant form of this ratio. We note that the PHOBOS observation [29] 
of a triangular dependence of the elliptic fiow V2{y) peaked at mid rapidity is very similar 
to the triangular pattern of {dETaii/ drj) / [dNch/ drj) predicted by HIJING due to mini-jet 
production at 130 A GeV. We are not aware of any predictions for this important global 
observable from saturation models. Fig. 3d suggests that the initial conditions for hydrody- 
namics are not well approximated by Bjorken boost invariant forms [41] assumed thus far. 
A full 3-l-lD hydrodynamical simulation [97] with such more realistic boost variant initial 
conditions should be investigated to try to account for the PHOBOS elliptic fiow. 

The PHENIX data [4] show a value closer to 0.8 GeV for {dExaii/di])/ {dNch/dr]) that 
is remarkably independent of ^sjvat from 17 GeV to 130 GeV and also independent of 
centrality. The observed independence on energy and centralities is very interesting since 
it is difficult to obtain such an effect in any transport theory with pQCD relaxation rates 
[69,100]. 

We investigate next in more detail the centrality dependence of dNch/dr) (Fig. 4) and 
dEx/dr] (Fig. 5) per pair of participating nucleons. Figure 4 and Fig. 5 show the results 
for centrahty dependence within HIJING vl.37 model calculations at y^Sivjv =130 GeV 
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and -y/sNN =200 GeV for (yq,ys) and (nq,ns) scenarios in comparison with experimental 
data [5,9,28,33]. The parameters employed for these calculations are shown on the figures. 
We note that all HIJING curves extrapolate at low multiplicities to the value dNch/dr]=2.2, 
observed in pp collisions by the UA5 collaboration [101]. The HIJING model predicts steady 
rise in the particle production per participant pair although the data seem to have a slower 
variation with Npart- The predicted increase is due to nonlinear increase of hard scatterings, 
which in contrast to the beam jet fragments, dependent on the number of binary collisions. 
HIJING/BB predicts a similar trend although the calculated values are lower than that 
given by HIJING by 10-15 % and under-predicts the experimental results at ^/sNN=130 
GeV. 

Figures 4a, 5a show the effect of lowering the energy loss to 0.5 GeV/fm as compared to 
a value of 2 GeV/fm used in Figs. 4b, 5b. The default parameter predictions at y^s/vjv=130 
GeV (Figs. 4b, 5b) are more consistent with data. However, for energy loss dE/dx =2.0 
GeV/fm and po— 2.0 GeV/c assumed to be independent of ^/s it is found that the ratio of 
i?2oo/i30 for midrapidity dNch/drj is over-predicted for most central collisions by 30% as shown 
in ref. [9]. This is a major failing of the HIJING assumption of energy independent mini-jet 
scale Po — 2 GeV/c. Motivated by the energy dependence predicted by the saturation scales 
in FSS and ISS models discussed in the introduction, and the data, we study in Figs. 4d and 
5d the effect of allowing a shght increase with energy horn po{y/s)— 2.0 GeV/c at y^sj^=130 
GeV to Po{a/s)— 2.18 GeV/c at ^/sNN— 200 GeV. Such an energy dependence was also found 
necessary in [54] using more modern structure functions than in HIJING. 

Figure 6 [Et/Nch transverse energy per charged particle) and Fig. 7 (ratios -R200/130 for 
midrapidity dNch/dr] and dET/drj) present the results obtained within both models HIJING 
vl.37 (upper part) and HlJlNG/BB vl.lO (lower part). The data are much closer to the 
quench and shadowing (yq,ys) scenario. 

It was shown [33] that hard scattering component to the charged particle production 
remains almost constant (36 ± 6) % over the energy range ^/snn=130 GeV - 200 GeV. We 
see when comparing the values from Figs. 4d,5d to Figs. 4c,5c and especially from Fig. 
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6 and Fig. 7, that the energy dependence of this global measure is reduced considerably 
to within the experiment range by allowing a modest increase of po without assuming any 
additional Npart dependence of this scale. From these results we conclude that a 10% in- 
crease with energy of the mini-jet scale (po) is required in both HIJING models to account 
for the centrality and energy dependence of the global multiplicity and transverse energy 
observables. 

III. JET QUENCHING AND THE NUCLEAR MODIFICATION FACTOR 

High pt hadron spectra have been widely analyzed at SPS and RHIC energies [40]. 
We investigate in this section how well can HIJING and HIJING/BB describe the high pt 
hadronic spectra in pp collisions and their predicted nuclear modifications in AA collisions. 

The observation [11,13,17,20] of strong suppression of high pt hadron spectra in cen- 
tral Au + Au at RHIC energies [79,81,82] is the most dramatic new dynamical phenomena 
discovered at RHIC relative to SPS. We recall that the comparison of parton model cal- 
culations and the experimental data does not show any evidence of parton energy loss at 
SPS energies [84]. The observed absence of quenching va. d + Au, [13,23,31,36] as predicted 
in Refs. [81,50,51,55,56], proves that quenching is caused by final state interactions in the 
dense matter formed in Au-|-Au collisions and not due to gluon shadowing. 

Parton-parton hard scattering with large momentum transfer produces high momentum 
quarks or gluons which fragment into jets of hadrons. The leading particles manifest them- 
selves in a power-law like shape of the momentum distribution. High momentum partons 
are predicted to lose a significant fraction of their energy by gluon bremsstrahlung leading to 
a suppression of the high momentum tail of the single hadron inclusive spectra [81]. It has 
been argued that data from RHIC experiments show characteristic features consistent with 
such "jet quenching effects" [6,105,106]. Other interpretations have been proposed after the 
data became available. These are based on gluon saturation in the initial nuclear wave- 
function [107], coherent fields and their geometry [108], surface emission of the quenched 
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jets [109], final state hadronic interactions [110] and quark coalescence [111]. 

The default HIJING implementation of jet quenching uses a simplified algorithm most 
closely resembling surface emission. The energy loss is implemented by testing the number 
of interactions that a jet will have along its propagation line with excited participant strings. 
The approximate linear participant number scaling of the bulk multiplicity motivates this 
approximation to the transverse matter density profile through which the jets propagate. 
The distance between collisions is fixed by a mean free path parameter, A = 1 fm by default. 
Energy loss is implemented by splitting the energy of the jet among multiple gluons with 
energies uji = AzidE/dx where dE/dx = 2(1) GeV/fm for gluon (quark) jets and Azi are 
distances between collisions. This simplified mechanism suppresses jets that originate more 
than one mean free path from the surface. 

The effect of the nuclear modification is quantified in terms of the ratio [53] 

<N^u>dmjdydp^' (1) 

where < Ncoii > is the average number of binary collisions of the event sample and can be 
calculated from the nuclear overlap integral (Taa) and the inelastic nucleon-nucleon cross 
section: < Ncou >— cr'^nn^ < Taa >■ In the HIJING effective surface emission model, we can 
expect Raa oc {X/Npi% fm). 

The absolutely normalized transverse momentum spectra and pseudo rapidity distribu- 
tions for Au+Au central (0-5%) coUisions [6,16] at ^/snn — 130 GeV are shown in Fig. 8. We 
compare the STAR [16] pseudo-rapidity distribution of negative hadrons (/i~) for the cen- 
tral (0-5%) Au-|-Au collisions with the predictions from HIJING (yq,ys)-solid line, HIJING 
(nq,ns)-dajsh.ed line and RQMD (dash-dotted hne). The negative hadron pseudo rapidity 
distributions in part (a) are best reproduced with shadowing and quenching effects present. 
Note in Fig. 8c, however, that the moderate 1 < pt < 4 GeV spectra are too strongly 
over-quenched by the default HIJING parameters. The "no shadow, no quench" calcula- 
tion over-predicts the global rapidity density but remarkably fits the moderate pt spectrum 
rather well. So where is the energy loss? As we discuss below, it is likely camouflaged by 
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anomalous baryon excess. 

Note that RQMD fits the moderate pt h~ data better than HIJING (yq,ys). However, as 
shown in Fig. 8b, the absence of mini-jet production in RQMD causes it to miss the observed 
features of the pp rapidity density [101] at -^snn— 200 GeV that HIJING reproduces. 
Similarly, because of multiple jet production, HIJING reproduces the power law like tail of 
the p+p pt spectra very well while RQMD shows a glaring discrepancy (Fig. 8d). This result 
is significant because it demonstrates that the agreement of RQMD in Fig. 8c is fortuitously 
due to a strong nuclear dependent Cronin like multiple collisions algorithm enhancement. 
This is similar to the fortuitous agreement [83] of HIJING with WA98 moderate pt pion 
data at the SPS, where the results were shown to be exponentially sensitive to the Cronin 
algorithm adopted. 

To better understand why HIJING over-predicts the quenching of h~ for < 4 GeV/c, 
we turn next to the latest data on identified 7r° central (0-10 %) Au-|-Au coUisions at -^snn— 
200 GeV. We use these data because of much higher pt reach than at i/sj^=130 GeV. Also 
in this case, Au + Au can be compared directly to the new p-|-p data measured by PHENIX 
[11,12]. In Fig. 9a, the predicted tt'^ pt spectra based on dE/dx—0, 0.5, and 2 GeV/fm are 
compared to the Au-|-Au data out to 8 GeV/c. In Fig. 9b the recent p+p — > n^+X data are 
compared to default HIJING. This shows even more clearly than the original test [37], how 
well HIJING is able to reproduce the high px spectrum out to 10 GeV/c in the elementary 
p + p case. The Raa{pt) nuclear modification factor is shown in Fig. 9c. It is observed that 
the default energy loss parameters(dE/dx=2 GeV/fm) describe very well the jet quenching 
pattern of neutral pions in central collisions. The Lund string fragmentation mechanism of 
hadronization in HIJING leads to a rather slow growth of Raa to unity at high pr (dashed 
histogram in Fig. 9c). Only after > 4 GeV/c do the details of hadronization become 
irrelevant, and in that range, the default energy loss leads to a factor of five suppression, in 
agreement with the naive surface emission estimate with A = 1 fm. The relative suppression 
ratio R2 = RAAiysq) / RAAinsq) (Fig. 9d) shows in more detail the levehng off of R2 beyond 
4 GeV/c for the default energy loss, while it levels off at 0.4 for a reduced energy loss of 
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0.5 GeV/fm. Note that in all calculations the default po — 2 GeV/c was used, as it has, 
however, no effect at high pt- 

In Fig. 10 we compare the quenching pattern of HIJING to data for peripheral (60-80%) 
reactions. Fig. 10b is the pp data scaled by the number of coUisions in this peripheral 
reaction class. The main difference with respect to Fig. 9 is that unhke in central collisions, 
even a reduced energy loss dE/dx = 0.5 GeV/fm over-predicts the small modification of 
Raa from unity observed for more peripheral interactions. As more clearly seen in Fig. lOd, 
HIJING predicts a 30% suppression in peripheral interactions that is not seen in the data. 
We conclude that while the central reaction suppression is correctly predicted by HIJING, 
the surface emission algorithm adopted to model energy loss is not realistic and does not 
reproduce the centrality dependence observed by PHENIX. 

We have studied the dependence of the quenching pattern on variations of the mean 
free path parameter (A) of HIJING as well. We find that the quenching is indeed sensitive 
to A. In central collision increasing A to 3 fm for example decreases the quenching for 
dE/dx = 2 GeV/fm by a factor of approximately two. To explain the centrality dependence 
an Npart dependence of A must be assumed. We do not pursue here such an elaboration of 
the HIJING energy loss algorithm but future studies would be desirable to help distinguish 
between surface emission and volume emission models of jet energy loss. 

The sensitivity of the nuclear modification factor Raa{p±) to the mean energy loss pa- 
rameter in HIJING at the lower 130A GeV energy is more clearly revealed in the high 
statistics computation shown in Fig. 11. The data at ^/sNN= 130 GeV are taken from 
PHENIX [6] and STAR [17]. This figure extends the comparison in Fig. 8 to the range 
Pt ~ 7 GeV/c. The main discrepancy between the negative hadrons h" data and calcu- 
lations is the observed distinct localized bump with a maximum at ~ 2 GeV/c. Raa 
approaches the predicted quenching pattern with dE/dx — 0.5 GeV/fm only at the highest 
Pt measured. 

PHENIX has found that the excess negative hadrons in the 2-4 GeV range are in fact 
due to anti-protons [3]. In order to check whether an enhanced baryon junction loop [86] 
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mechanism could possibly account for that excess, we plot Raa for HIJING/BB [87] in Fig. 
12. While the junction source reduces the discrepancy between data and HIJING shown 
in Fig. 12a, it fails to account for the very large excess of anti-protons at moderate pr- 
From Fig. 12a,b, we see that junctions as currently implemented do not solve the centrality 
dependence problem discussed in connection with Fig. 10 either. 

As a final comment we compare the "central to peripheral" nuclear modification factor 
defined by 

pep / N _ {Yield/ < N^ou >)(o-5%) 
""^^^^^^ {Yield/ < N^oii >)(60-80%) 

where Yield — {l/Nevents){^/'^'^P±){d'^N/dp±dr}) to calculated values for this particular ra- 
tio. The data are from PHENIX [6,8] and STAR [17]. Even though BB version fails to 
describe both the numerator (Fig. 12a) and the denominator (Fig. 12b), it accidentally 
describes the it^Xi ratio of central collisions (Fig. 13a). No such lucky coincidence occurs 
for peripheral coUisions using the default HIJING (Fig. 13b). This figure demonstrates 
the great care one must exercise in interpreting any agreement of dynamical models with 
specific ratios. It is always essential to check whether the model is able to reproduce the 
absolutely normalized spectra, as in Figs. 9 and 10. Only after a model passes that test can 
any agreement with specific data ratios be considered seriously. 



IV. SUMMARY AND CONCLUSIONS 

We have investigated in this paper how the predictions of HIJING and RQMD exclu- 
sive nuclear collision event generators compare to the new available data from RHIC. We 
concentrated on two classes of observables. First the global number and transverse energy 
distribution in rapidity was considered. Then we focused on the new jet quenching nuclear 
modification factors. 

The energy dependence of global observables rule out RQMD because of its neglect of 
hard pQCD mini-jet production. However, the observed energy dependence also rules out 
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HIJING in its default parameter settings. The separation scale, po between soft and hard 
processes, which assumed in HIJING to be a constant 2 GeV/c independent of ^/s and 
centrality, predicts a too rapid growth of multiplicity. Motivated by FSS and ISS parton 
saturation models and the data, we tested and found that allowing a 10% growth of po from 
2.0 to 2.18 GeV greatly improved the consistency of HIJING results with the observed RHIC 
systematics. In all cases, the default shadowing (with identical quark and gluon shadowing) 
assumed was found to be essential to reduce the mini-jet flux and not to over-predict the 
multiplicity. The small enhancement of the multiplicity due to jet quenching with the default 
energy loss was consistent with experimental data once the small energy dependence of po 
is taken into account. 

Our analysis of the jet quenching pattern predicted by HIJING shows that the default 
dE/dx = 2 GeV/fm accounts remarkably well for the suppression pattern of tt*^ out topr = 8 
GeV/c as observed for central Au + Au coUisions at y^sjviv=200 GeV. A major advantage 
of HIJING over other models is that it reproduces accurately both the low pt dominated 
rapidity density and the high pt recent 200 GeV ■n'^ + X data [12] at the same time. 

However, neither HIJING nor HIJING/BB are able to account for the anomalous baryon 
lump in the intermediate < 4 GeV/c region. Furthermore, we noted that the energy 
loss algorithm in HIJING corresponds effectively to surface emission with a default A = 1 
fm mean free path. We checked that increasing A leads to less suppression. We found that 
a constant N^art independent A, however, is not compatible with the observed centrality 
dependence of jet quenching. 

The failure of the current implementation of baryon junction loops in HIJING/BBvl.lO 
to reproduce the observed pr enhancement of anti-baryons and baryons needs, however, 
further study since an enhancement was theoretically anticipated [89,86]. We are currently 
investigating why this feature of baryon junction dynamics did not emerge from numerical 
simulations with this code. Understanding the physical origin of the (anti) baryon anomalies 
is essential to disentangle competing mechanisms such as collective hydrodynamic flow [41], 
multi-quark coalescence [111], and possibly novel baryon junction dynamics [48] at RHIC. 
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FIG. 1. Charged hadron rapidity distributions for central (0-6%) Au+Au collisions as function 
of cm. energies. The histograms show the theoretical predictions from HIJING vl.37 (upper 
part) and HIJING/BB vl.lO (lower part) with (y) or without (n) effects of quenching (q) or/and 
shadowing (s) included. The data are from PHOBOS collaboration [26] (a), [27] (b) , [28] (c) and 
BRAHMS Collaboration [33] (c). The error bars at midrapidity include systematic uncertainties. 
The others error bars of the order of 10-15 % have been omitted for clarity. 
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FIG. 2. Charged particle rapidity density per participating baryon pair versus the cm. energy. 
The predictions of HIJING/BB vl.lO (yq,ys)-Ml line, HIJING/BB vl.lO fn^^ns^-dashed hne, HI- 
JING vl.37 (yq,ys)-dotted line, HIJING vl.37 (nq,ns)-dot dashed line, and RQMD v2.4 {open 
squares) are compared to data. The data for the central (0 — 6%) Au+Au collisions, are from 
PHOBOS [26], [28], and from BRAHMS [33]; the error bars include systematic uncertainties, pp 
and pp data are from ref. [101-104]; the error bars are statistical only. 
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FIG. 3. Transverse energy distribution for: (a) all particles, (b) for charged particles only, (c) 
total transverse energy per charged particles and (d) total transverse energy of charged particles 
per charged particles as a function of pseudo rapidity. Theoretical predictions from HIJING/B B 
vl.lO (yq,ys)-soM, HIJING vl.37 (yq,ys)-dashed, HIJING vl.37 (nq,ns)-dash dotted and RQMD 
v2.4-dotted histograms are shown. The data are from PHENIX [5] and STAR [16]. The error bars 
for PHENIX include systematic uncertainties. 
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FIG. 4. Midrapidity dNch/drj per participant pair as a function of the number of participants at 
y^SiViV=130 GeV (upper part). Theoretical predictions from HIJING vl.37 model with (ysq-solid 
lines) and without (nqs-dashed lines) the effects of quenching and shadowing. The data at 
y's]^= 130 GeV are from PHENIX [5], [9], PHOBOS [28]. Lower part are the results at ^s]^=200 
GeV. The data are from BRAHMS [33], PHOBOS [28] and PHENIX [9]. The parameters within 
HIJING calculations are given in the figure. The error bars include both statistical and systematic 
uncertainties. 
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FIG. 5. Midrapidity dEx/dr] per participant pair as a function of the number of participants 
at y^sivjv=130 GeV (upper part) and y^sjvjv=200 GeV (lower part). Theoretical predictions from 
HIJING vl.37 model with (ysq-solid lines) and without (nqs-dashed lines) the effects of quench- 
ing and shadowing. The solid and dashed lines have the same meaning as in Fig. 4. The data 
at ./sWn=130 GeV are from PHENIX [5], [9], PHOBOS [28] and at ^/sWn=200 GeV are from 
BRAHMS [33], PHOBOS [28] and PHENIX [9]. The error bars include both statistical and sys- 
tematic uncertainties. 
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FIG. 6. Midrapidity ratios {{dET / difj / dNch / drj) as a function of the number of participants at 
V*iVAr=130 GeV and y^sivjv=200 GeV Theoretical predictions within HIJING vl.37 (upper part) 
and HIJING/BB vl.lO (lower part). The solid and dashed lines have the same meaning as in 
Fig. 4. The data are from PHENIX [9]. The error bars include both statistical and systematic 
uncertainties. 
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FIG. 7. Ratios -R200/130 for midrapidity dNch/drj and dET/drj per participant pair as a function 
of the number of participants. Theoretical predictions within HIJING vl.37 (upper part) and 
HIJING/BB vl.lO (lower part) model with (ysq,po) a-nd without (nsq,po) the effects of quenching 
and shadowing are compared to data [9]. For the labels of solid and dashed lines see Fig. 4(b), (d) 
and Fig. 5(b),(d). 
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FIG. 8. Comparison of Au+Au and p+p at RHIC energies is shown, (a): Negative hadrons 
pseudo-rapidity spectra for Au+Au collisions at y^sjvjv=130 GeV including STAR central (0-5%) 
data [16]. (b) p+p at y^s]^=200 GeV including UAl p+p data [101]. HIJING vl.37 (yq,ys)-soM, 
HIJING vl.37 (nq,ns)-dash.ed, and RQMD v2.4-dot dashed histograms are theoretical calculations. 
Parts c) and d) compare the p± distributions. PHENIX data [6], STAR data [16] and UAl data 
[104] are shown. The error bars include systematic uncertainties in part (a) and are statistical only 
in parts (b,c,d). 
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FIG. 9. Invariant yields at midrapidity for 7r° in central (0-10 %) Au+Au collisions (Part a) 
and scaled by number of binary collisions (Ncou) in p+p interactions at ^/sNN =200 GeV. Nuclear 
modification factor Raa (Part c) and ratio R2 (Part d) as a function of transverse momentum 
as predicted by HIJING vl.37 with (solid and dotted histogram) and without (dashed histogram) 
shadowing and quenching effects. The labels have the same meaning as in Fig. 4. The data are 
from PHENIX [10-12]. Only statistical error bars are shown. 
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FIG. 10. Invariant yields at midrapidity for 7r° in central (60-80 %) Au+Au collisions (Part a) 
and scaled by number of binary collisions (J^coll) in P+P interactions at y^SiV7V=200 GeV. Nuclear 
modification factor Raa (Part c) and ratio Ki (Part d) as a function of transverse momentum as 
predicted by HIJING vl.37 with (solid histogram) and without (dashed histogram) shadowing and 
quenching effects. The labels have the same meaning as in Fig. 4. The data are from PHENIX 
[10-12]. Only statistical error bars are shown. 
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FIG. 11. Nuclear modification factor Raa, as a function of transverse momentum as predicted by 
HIJING vl.37 for central (0-5%) Au+Au collisions at ./snn=130 GeV . The different histograms 
are the results obtained without quenching and with quenching assuming different values for parton 
energy loss. The data are from PHENIX [6] (squares) and from STAR [17] (stars). Only statistical 
error bars are shown. 
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FIG. 12. Comparison of the nuclear modification factor (Raa) for Au+Au as predicted by 
HIJING/BB vl.lO (upper part) and HIJING vl.37 (lower part) for central (0-5%) and peripheral 
(60-80%) Au-FAu collisions. The data are from PHENIX [6,8] and STAR [17]. Only statistical 
error bars are shown. 
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FIG. 13. Comparison of the nuclear modification factor (-RX4) Au+Au as predicted by 
HIJING/BB vl.lO (part (a)) and HIJING vl.37 (part (b)) for ratio of central (0-5%) to peripheral 

(60-80%) collisions. The data are from PHENIX [6,8] and STAR [17]. Only statistical error bars 
are shown. 
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